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ABSTRACT 

The Amazon Region natural resources are 
studied in two ways and compared. Initially, a 
LANDSAT scene and its attributes are selected, 

* ami a maximum likelihood algorithm classification 
Is made. The next step is the atmospheric 
correction of the scene, taking into account 
Amazonic peculiarities (checked with ground 
truth) of the same area, and the subsequent 
classification. Finally a comparison, wich 
shows {hat the classification improves with the 

. atmospherically corrected images, is made. 

# * 1 INTRODUCTION 

f In order to better use orbital remote 

' sensing techniques for natural resources studies, 

' it is desirable to minimize the nontarget-related 
3 attributes influence, as that of atmospheric 
! absorption and scattering. 

* Therefore, the main objective of this 

: paper is to analyse the atmospheric influence on 
LANDSAT imagery for the Amazon region data. 

i 

{ The selected area was centered around 

the city of Manaus (Capital of Amazonas State, 
Brazil), containing the Manaus* International 
t Airport (Eduardo Danes Airport), from which 
radiosonde data are available. 

The results will be used at INRE 
(Institute for Space Research, Brazil) and other 
Brazilian institutions that utilize remotely 
sensed LANDSAT data for Amazon Region studies. 

II - STUDY AREA 

The selected study area is the scene 
sub-image 248.62 (WRS). The land area is 

* approximately 900 km*, located between 

$V> 50* W - 60° 00* W and 03° 00* S - 03° 10* S. 
Within the selected area the city of MANAUS, its 
’ INTERNATIONAL AIRPORT, jungle plantations, and 


rivars with different types of water, which 
makes it a good choice, are located. The proximity 
of MANAUS enable one of the authors to collect 
data >n site for use as ground truth. 

The LANDSAT material used were CCT*s in 
four channels, for September 28th, 1981, and the 
machine processing performed on INPE*s IMAGE-100 
System,, assembled by GENERAL ELECTRIC CO. 

• III - METHODOLOGY 

A procedure developed at INPE for 
atmosoheric correction 1 was used. Its steps are 
* as follows: (i) the atmospheric transmittance is 
determined by a modified LOWTRAN 4 2 computer 
program (originally developed by the Air Force 
Geophysical Laboratory, AFGL, USA). This program’ 
was modified for interactive remote video 
terminal use, with optional graphical output*. 

The atmospheric model, inferred from radiosonde 
data, was acquired at Eduardo Gomes Airport 
(Manaus, AM, Brazil), at the same date and 
approximately at the same time of the satellite 
pass; (ii) using the transmittance, the solar 
zenith angle, and the estimated values for the 
average ground albedo and aerosol optical depth, 
the LANDSAT image is corrected pixel by pixel 
by the 1-100 System. 

♦ In order to analyse quantitatively the 

atmospheric correction effect, it was necessary 
to obtain true grey level range for the selected 
classes identified in the scene. 

The procedure,, with and without 
atmospheric correction, for obtaining the 
spectral classes was the following: (i) noise 
removal, (ii) radiometric correction to minimize 
stripping, (iii) main area selection at 1:100,000 
scale, (iv) training areas selection (obtained 
by field trip ground truth); (v) spectral 
parameters acquisition through the MAXVER 4 
(maximum likelihood) algorithm. 
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Other quantitative parameters were 
measured in order to compare the data with and 
without atmospheric correction. Initially* the JM 
distance (Jeffreys Matusita Distance) was computed 
to determine the separability among the spectral 
classes for one or a combination of channels. The 
JM distance 5 is given by the formula: 

JM - 2 (1 - e“°) , 
where 4 

or ■ i (Ui-Uj) e * 1 (U!-U 2 )*i In 

e 2 

Cv 1 ... i. ; it L v- ..v.-i * 


(D 


jgj 6 

/ det ci det 


and 


U, ■ averages of class 1 matrix 0, 1 • 1,2;* 
• * T ‘ 

( ) * transposed matrix; 

( )"* • inverse matrix 

e • - [cj ♦ ei]; 

2 

c, ■ covariance matrix class 1, 1 ■ 1,2; 

*1 

dete« determinant, 1 * 1,2. 

The applicability of these measurements 

for class separation, mainly for vegetatiorf 
’•covered areas, was shown by Aoki and Santos 6 * 

Another useful parameter is the ratio 
between the average grey levels with and without 
• atmospheric correction, for each class, for each 
MSS channel. This ratio is expressed by: 

ra(i, 6) . ills®? , O) 

f C(1 ,6) 

being C'(i,B) > C(i,S), 

where C'(i,B) * average gray level for class 1, 
channel 6, with atmospheric 
correction. • • 

C(i,8) » average gray level for class i, 
channel S, without atmospheric 
correction. 

Next, the product PA, composed of the 
contrast index between two classes multiplied by 
the average value of these classes is studied: 

.. £!idl , ( «) 

• C(J,9) 2 

. t c-xrr-r.ii t;< ;*<•••»• • :i t s ' :no I 


f 






for 1 > j, 

being C(i,6)/C(j,6) the contrast Index, and 
(C(1,B) c(j,3})/2 the average, where C(i,8) ■ 

average gray level value for class 1 (or j), 
without atmospheric correction; and 8 is LANDSAT 
MSS channel. 

The product PA 1 Is similarly obtained 
for atmospheric corrected values. 



Following, the ratio RB (1,j,B) between 
PA* and PA for PA* > PA is: 

RBOJ.B) * PA'(i ,j ,8)/PA(1 ,j ,B) . (5) 

Finally, the standard deviation for 
each class, for each MSS channel, with and 
' without atmospheric correction, taking into 
account the size of the training area, is 
. analysed. 


IV - RESULTS 

Fran radiosonde data, assuming a 
visibility of 20 km, the atmospheric transmission 
was computed through the use of a LOWTRAN 4*: 

Ti, - 0.6905, 

T* - 0.7467, 

T,.« 0.7343, 

T, ■ 0.7646. 

The zenith angle (50°), aerosol optical 
depth assumed (0.12), and estimated albedo (0.20) 
are the other input parameters for the 
atmospheric correction program 1 . A composition of 
channels 4,5 and 7, atmospherically corrected 
is shown in Figure 1. 

It is noted on the picture that the 
Negro River water and the Solimoes River 
(Amazonas River) water are strikingly different. 
The city of Manaus is in the middle of the 
picture. This picture is a sub-image of orbit 
248.62 (WRS), with a 1:100.000 scale on the 
1-100 video. 



1983 Machine Processing of Remotely Sensed Data Symposium 


noiRlN&L PAGE {9 



Figure 1 - Area under study with atmospheric 
correction. 

In the selected scene, six classes were 
considered: (1) dark water, (2) muddy water, 

(3) jungle (tropical rain forest), (4) secondary 
forest, (5) urban areas (city of Manaus and 
. environs), and (6) agricultural areas. The 
classes spectral parameters with and without 
atmospheric correction, obtained by the MAXVER 
algorithm, are presented on Table 1. 

The values of RA for each class and 
channel are presented on Table 2, 

Table 3 presents RB for each -pair of 
■ classes, while Table 4 presents the standard 
deviation as computed by MAXVER. 

V - DISCUSSIONS AND CONCLUSIONS 

From Table 1, it is observed that, 
depending on the target, the correction in 
general increased the average grey levels. For 
example for class 1, cnannel 7, the difference 
was 118*, while for class 4, channel 4, the 
difference was only 25 ?. The value of the 
correction is highly dependant on the target. 

It is observed from Table 2 that class 1 
was more affected by the correction, followed by 
class 2, while on the other classes the effect 
was not so prominent, 

Tabl." 3 shows that the RB value did not 
change very much for pairs of classes on the 
visible and mar infrared (channels 4,5,6), while 
for the infrared channel 7 it varied from 1,00 
to 1,66. 


The standard deviation increase with 
the atmospheric correction; thus, it is supposed 
that a class has a larger set after correction, 
which should improve the classification. It 
! should be noted that class 3 (jungle) presented a 

■ much larger standard deviation than the other 
: classes, while class 4 {secondary forest) 

presented a small standard deviation. This result 
J is still under study. 

In addition to these quantitative 
results, it was observed that after an 
' atmospheric correction is performed there is a 
noticeable increase in the image visual quality, 
with sharper natural ground target contours. 
Finally, the above results are preliminary, and 

■ many more tests and field trips are needed to 
achieve a better understanding of the atmospheric 

■ effect on LANDSAT imagery over the Ama2on Region. 
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